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ABSTRACT 


Warping constraints of thin-walled beams of open cross sections may 
significantly add to the torsional rigidity of the beams. This property is 
not included in the conventional beam element available in NASTRAN and makes 
it impossible to model such beams properly. 

This paper presents a composite element which includes the effects of 
warping and offset shear centers. It is shown to be mathematically consistent 
with thin-walled open beam theory and can easily be incorporated into any 
structural analysis program by use of NASTRAN' S standard elements. 

A numerical example analyzing the vibration of a channel beam is presented, 
and, by using the composite element, the numerical results agreed very well 
with theoretical data. 


INTRODUCTION 


Thin-walled beams of open cross sections are commonly used in automotive 
and aircraft structures where it is imperative to insure maximum flexural 
efficiency without violating practical fabrication and assembly requirements. 

It turns out that for many cases the static and dynamic response of these 
structures depends strongly upon the warping constraints of the members. Con- 
sequently, proper modeling techniques are essential to insure good numerical 
results . 

The most general beam element available in the present version of NASTRAN 
(level 12) is based upon the well-known Timoshenko beam theory. This theory, 
however, neglects warping displacements and assumes that the elastic and 
centroidal axes coincide. Both these conditions are usually too restrictive 
when thin-walled open beams are modeled. 

Several books and articles have been published treating the theory and 
applications of thin-walled beams of open cross sections. A thorough presenta- 
tion is found in Vlasov's book (Ref. 1). Timoshenko and Gere have done a 



buckling analysis (Ref. 2), and contributions to the vibrational analysis 
are due to Timoshenko (Ref. 3), Gere (Ref. 4), Gere and Lin (Ref. 5), and 
Chris tiano and Salmela (Ref. 6). 

This paper presents a new element which includes the effects of warping 
and noncoincidal elastic and centroidal axes. It is a composite element 
derived from a variational principle consistent with the thin-walled open beam 
formulation. Also, in addition to being mathematically correct, this model can 
very easily be incorporated into any structural analysis program by using 
elements now available in NASTRAN. 

As a numerical example on how the model is constructed, the vibration of 
a channel beam is analyzed, and good agreement between theoretical and numeri- 
cal results are obtained. 

THEORY OF THIN-WALLED OPEN BEAMS 


A thin-walled beam is usually defined by restricting its dimensions so 

that 


L > 10d 

(1) 

d > lOt 

(2) 


where L is the length of the bar, d the depth of the cross section, and t the 
largest thickness (Fig. 1). The first relation makes it possible to assume 
uniformity of the longitudinal stresses along the axis (St. Vernant’s 
Principle) while Eq. 2 is used in making assumptions on the shear stresses in 
the section. Particularly for beams of open cross sections which will be 
considered in this analysis, it will be assumed that the middle surface is 
free of shear. In addition, the material and sectional properties of the beam 
are assumed to be constant along the length, and the displacements are con- 
sidered to be small such that the cross sections do not change shapes during 
deformation. 

Figure 1 shows a general thin-walled open beam. Included is also a 
detail of the cross section giving the locations of the center of gravity (CG) 
and the shear center (SC) . The shear center is the point on the cross section 
through which the resultant of the transverse shearing forces always passes. 
Locus of the shear centers along a beam is called the elastic axis, and 
similarly, the centroidal axis passes through the CG of the cross sections. 

Before the differential equations can be presented, it is very important 
to clearly define the coordinate and displacement system used. Referring to 
Fig. 1, the x-coordinate coincides with the centroidal axis of the beam, and 
the y- and z-coordinates are the principal axes of the cross section. The 
displacements are defined by 

u « longitudinal displacement of centroid 

v = vertical displacement of shear center 
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w = horizontal displacement of shear center 
0 = rotational displacement of shear center 

and are also indicated in the figure. 

Based on these definitions and the assumptions mentioned above, the 
full set of differential equations for a thin-walled open beam can be written 
as (Ref. 1) 


•pa d 2 u 
EA - q 

dx 2 x 


(3) 


El! ^ - q 
dx 4 y 


(4) 


EI 2 ^ = q 
dx 4 Z 


(5) 


EC ££ _ G J — = m 
W dx 4 dx 2 9 


( 6 ) 


in which E = modulus of elasticity; G = shear modulus; A = cross sectional 
area; Ii = moment of inertia about z axis; I 2 = moment of inertia about y-axis; 
J = torsion constant; Coj = warping constant; q x = longitudinal surface load 
applied to the centroid; qy and q z = transverse loads applied at the shear 
center; and mg = twisting moment about the elastic axis. 

It should be observed that the above set of differential equations 
differs from the regular beam equation by the fact that v, w, and 0 (and 
correspondingly q y , q z , and m Q ) are displacements (and forces) relating to 
the shear center rather than the centroid , and that the term 


dx 4 

relating to the warping of the cross section, has been added. Also, the 
four equations are completely uncoupled, a very desirable feature made pos- 
sible by letting some of the displacements be associated with the shear center. 

MODELING OF THIN-WALLED OPEN BEAMS 

The conventional beam element available in the NASTRAN computer program 
is essentially based upon the same assumptions given above, but in addition, 

two other conditions are made: 

1. Plane sections remain plane (i.e. , no warping deformations) 

2, Centroidal- and elastic-axes coincide. 
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These two additional restrictions were also pointed out in the last section 
by looking at the differential equations . 

Warping is a phenomenon that occurs in all bars subjected to twist, 
with the exception of circular cross sections (Ref. 7). Only for open 
sections though can the strain associated with the warping be large enough 
to cause significant contribution to the rotational stiffness of the beam. 

By excluding warping, the rotational stiffness is too flexible; a factor of 
two or more is quite common depending upon the cross-sectional properties 
and how the beam rotates. There have been some attempts to incorporate warp- 
ing in the conventional beam element by defining an effective torsional con- 
stant J e . From Eq. 6 it is observed that this can be accomplished by the 
equation 


GJ de _ de 

Gj e dx “ GJ dx 


EG 


d 3 e 

dx 3 


(7) 


but it is also noticed that J e then necessarily must depend upon the solution 
0. For static analysis, this approach might be used effectively where the 
rotational displacement 6 can be estimated fairly accurately, but for dynamic 
analysis, where different modes are extracted in the same analysis, large 
errors will be introduced. 

When centroidal and elastic axes do not coincide, the beam element 
should be modeled along the elastic axis. Otherwise, it can easily be observed 
that the bending moment created by the applied transverse forces will not be 
encountered. This modeling will unfortunately introduce new problems, namely, 
the bending moment caused by the longitudinal forces which, according to Eq. 3, 
should coincide with the centroidal axis. In many cases, though, the effect 
of axial forces can be neglected. 

In order to treat thin-walled open beams in a finite element computer 
program and include the effects of warping and offset shear center, it will 
be necessary to either 1) develop a new element that includes these effects, 

2) model the beams by a number of plate elements, or 3) arrange the avail- 
able finite elements in a way consistent with the thin-walled open beam 
theory. While there exists an element (Ref. 8) with the above features, it 
is noc readily available, and in particular, not for NASTRAN users. A beam 
modeled by plate elements usually increases the number of grid points tremen- 
dously. This also increases the degrees of freedom, and the procedure will 
often become impractical and uneconomical. One is thus often left to model a 
thin-walled open beam by the existing beam element which will introduce errors 
as discussed above or create a new composite element* The following section 
presents such an element where the effects of warping and offset shear centers 
are included. 
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DERIVATION OF A COMPOSITE ELEMENT 


Before describing a suitable finite element model, a fictitious continu- 
ous beam is presented (Fig. 2). The differential equations of vibration 
and the corresponding boundary conditions for this beam will be derived 
by the Hamilton's Variational Principle 


6 


Jt 

I 


(T - U) dt - 0 


( 8 ) 


where T is the total kinetic and U the total potential energy of the system. 


Referring to Fig. 2, the fictitious element consists of a rod coinciding 
with the centroidal axis, a beam coinciding with the elastic axis, and two 
flanges displaced h from the elastic axis in a direction parallel to the 
y coordinate. The rod accepts only axial deformations, and the corresponding 
strain energy over the whole length, L, of the rod can thus be written as 


U 


R 




(9) 


where u is the longitudinal displacement of the centroid as defined in the 
last section. 


As indicated in the figure, the beam has all capabilities of a regular 
beam except that the area is neglected such that no axial strain exists. 

The strain energy is then given by 



( 10 ) 


in which I 2 * is the moment of inertia about the y-axis while the other 
properties, including the displacements, are as defined earlier. 

The two flanges are displaced from the elastic axis an arbitrary dis- 
tance h which will be determined later. Only bending deflections are admitted 
for the flanges, and by defining the displacements parallel to the z-axis as wi 
and w 2 , as indicated on Fig. 2, the strain energy becomes 
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The rod is assumed to carry all the mass of the elements such that the 
kinetic energy will he written in terms of the rod displacements u, v*, w*, 
and 0*. Neglecting rotational inertia, the kinetic energy is 



(12) 


where Ip is the centroidal polar moment of inertia and p is the mass density 
of the material. 


The nine displacement parameters used will next be reduced to four by the 
relations 


8* 

= 0 



(13) 

V* 

= V 

+ 

c 0 
z 

(14) 

w* 

= w 

- 

c 0 

y 

(15) 

Wl 

= w 

-f 

h0 

(16) 

W2 

= w 

- 

he 

(17) 


Cy and c z are the distances between C.G. and S.C. in the y and z direction 
respectively. Substituting these equations into the above expressions for 
strain and kinetic energies, Hamilton's equation (8) can be formed. The 
variation is then accomplished by varying u by 6u, v by <5u, w by <$w and 6 by 
60 where it is understood that the variation vanishes at t = t 0 and t = tj. 
This procedure produces the differential equations 


3 2 u . 3 2 u 

EA — = pA — - 

3x 2 3t 2 


(18) 


3x 4 3t 2 


pAc 


3t 2 


(19) 


E (l 2 * + 2I F ) 


3 4 w _ ^ 3^ 3 2 0 

3x 4 ” P 3t 2 P ^ 3t 2 


( 20 ) 


2E%h 2 — - GJ — 
3x 2 3x 2 


,2 + c 2 + i - 


~ pA ic z A + c y ' a/ 2 


pAc_ 


3 z v 


P ACy 


3^ 


3t 2 7 3t 2 


( 21 ) 


and the boundary conditions (at x = 0 and x = 1) 
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( 22 ) 


afv 

. 3v . 

3x 2 

3w 

3 3 v 

• 

3x 3 

V = 

3 2 w 

• . 

Bx 2 

3x 

3 3 w 

3x 3 

. w = 

3u 

9x ’ 

U = 1 

3 2 e 

. 90_ ; 


3*2 3X 

(2EI_h 2 — - GJ —V 9 = 0 
V F 9x 3 W 


(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


The differential equations presented for the thin-walled open beam, 

Eqs. 3 through 6, can easily be transformed into the form given for the ficti- 
tious beam by substituting inertia loads for q x , q y , q 2 , and m Q (Ref. 1), 
Equations 18 through 21 will then become identical with the thin-walled open 
beam equations if the following two substitutions are made 


1.2* = Iz - 2 Ip 


(29) 


I 


F 



(30) 


It should be observed that the boundary conditions also will be identical 
with these substitutions. 


Because of the mathematical identity between the fictitious and thin- 
walled open beams, a finite element model for the first will also be a model 
for the latter. Figure 3 shows how the fictitious beam can be modeled and how 
the different elements are constrained at the grid points a D and b Q . Four 
elements will be used; a rod, without torsional rigidity, coinciding with the 
centroid; a beam, without axial rigidity and mass, coinciding with the elastic 
axis; and two flanges (or beams), resisting only bending deformations, and 
displaced h parallel to the y-axis away from the elastic axis. The second 
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moment of inertia (NASTRAN notation) for the beam and the flanges is as 
defined in Eqs. 29 and 30, otherwise the notation given in Fig. 3 refers to 
the original properties of the thin-walled beam. The distance h is arbitrarily 
chosen, but should he large enough such that I 2 * is positive. 

At each grid point, which is arbitrarily taken to be the centroid, the 
beam is offset by a rigid link to its shear center location. This will insure 
that the forces acting on the grid points will be transmitted correctly 
to the beam elements, and the grid point displacements will obey Eqs. 13, 14, 
and 15 used in the development of the fictitious beam. 

Similarly, the displacements of the two flanges must be made dependent 
upon those at the shear center as restricted by Eqs. 16 and 17. If warping 
is constrained at the grid point, i.e. , 


de 

dx 


0 


(31) 


then the connections between the flanges and the centroid can be made rigid. 
Otherwise, as in the case of interior points of a thin-walled open beam, 
warping is admissible, and Eqs. 16 and 17 imply that 


dx dx’ dx dx 


(32) 


This makes it necessary to create special grid points (aj, a 2 , bi and b 2 
on Fig. 3) for the flanges where the transverse displacements are made depen- 
dent displacements, and the rotational displacements 


dwi dw 2 

j — and j — 
dx . dx 

are unrestricted. The four other displacements must, in this case, be 
constrained in order to avoid singular matrices. 

The only difference between the composite model and the fictitious 
beam is that while Eqs. 16 and 17 are satisfied continuously for the mathe- 
matical beam, they are satisfied only at discrete intervals for the model. 
This implies that the flanges between the grid points will deform according 
to the law of minimum energy, contributing less to the torsional rigidity 
than theoretically assumed. Although usually small, this error can be made 
still smaller by using more grid points along the length of the beam. 


NUMERICAL EXAMPLE 


A numerical example is presented next to illustrate the use of the com- 
posite element and to compare the numerical results with data obtained from 
theory. The example chosen is the dynamic response of a channel beam whose 
properties are given by 
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A =0.884 in 2 p = 0.733 * IQ" 3 lbs sec 2 in- 4 

11 = 7.66 in 4 C = 3.52 in 6 

1 2 = 0.294 in 4 c“ = 0.0 in 

J = 0.00168 in 4 e y = 0.94 in 

E = 29 * 10 6 lbs in- 2 L Z = 120 in 

G = 11 • 10 6 lbs in" 2 

In the pure vibration of the channel beam, the differential equation in the 
axial displacement u (Eq. 18) uncouples from the other equations. Since Cy 
is zero in this example, Eq. 20 will uncouple also. The problem then reduces 
to a coupled vibrational problem in the v-* and 0- displacements, Eqs. 19 and 
21. For simply supported ends. 


v 



e 


i!i. o 

ax 2 


(33) 


and the theoretical solution, as given by the natural frequencies of vibration, 
p±, can be written as (Ref. 3) 


2 _ 


0) 


ti 


+ to. 


bi 




to . 2 — 

ti 


bi 


4 Ac to, 
z bi 


“ti 


2(1 - 


Xo z> 


where 



o 


(34) 


(35) 


I = I + A (c 2 + c 2 ) (36) 

o p \y z / 

and to t i and to^ are the i cn uncoupled frequencies in twist and bending, 
respectively; 



GJir 2 L 2 i 2 + EC ir 4 i 4 

to 

pi L 
o 


(37) 



Eliir 4 

i 4 

PAL 4 


When warping is not included in the analysis, to t i reduces to 


(38) 
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GJtt 1 2 3 L 2 1 2 


(39) 


The data obtained by using these equations are given in Table 1, cases 
1 through 4. It can be observed that warping increases the twisting frequen- 
cies more than 100% while an increase of about 5% is seen for the bending 
frequencies by going from uncoupled to the coupled theory. 

The channel is modeled by five equal elements 24 inches in length. Seven 
different modeling techniques have been tried, and the result is given in 
Table 1. Cases 5 through 8 show the channel modeled by the conventional 
Timoshenko beam element coinciding with the centroidal or elastic axis. 

In the first case only the uncoupled frequencies are obtained, while case 8 
reflects the coupled frequencies without taking warping into effect. Also 
noticeable are the bending frequencies obtained in case 7, which because of 
the coupling between v and 9, actually represents the twisting frequencies. 

Warping is included in the last three cases in Table 1 by using a com- 
posite element similar to the one given in Fig. 3, but the rod element is 
excluded since no axial forces are present during the vibration considered. 
Because of the boundary conditions, Eq. 33, the flanges must be free to 
rotate at the ends. This is accomplished by modeling the ends similar to 
all the interior points with a semi-rigid link where w^ and W 2 are restricted, 
but the rotations about the y- axis are free. The bulk data cards used for 
case 11 are given in the Appendix. Case 9, where the flanges are rigidly 
linked to the shear center, is included to show how important it is to free 
the rotational coordinates of the flanges, while cases 10 and 11 show the 
composite element properly used, including bending displacement and both 
bending (v) and twisting (8) displacements, in the numerical analysis. It is 
noticed that good correlation is obtained between theoretical (case 4) and 
numerical (case 11) data, the latter results being slightly below the first. 
This can be reasoned by the fact that few elements are used in the analysis, 
which implies 1) reduced torsional rigidity as explained above, and 2) con- 
centrated inertia forces at the grid points instead of a continuous 
distribution. 

Although the example did not use the rod element, this must be considered 
a very special case. For other structures and other loading conditions, 
the effect of excluding the rod can be significant. This will especially 
be the case where thin-walled open beams are used as members in between closed 
beams and where the axial load transfer is large. The example has shown, 
though , that 

(1) warping effects are properly taken into account by the 
composite element 

(2) inclusion of rotational displacements are necessary 

(3) good results are obtained with few elements. 
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CONCLUSION 


The composite element presented in this paper has been shown both mathe- 
matically and numerically to properly represent the thin-walled open beam. 

Since it is derived upon the differential equations for such beams, the model 
can be made as accurate as the theory is. In fact, in the limit as the length 
of the element gets smaller and smaller, the composite beam element is an exact 
model of the thin-walled open beam. 

In conclusion, it is hoped that future editions of NAS TRAN can present 
a true thin-walled beam element since the need for proper modeling will 
increase. Until that time, the present composite element is available, it 
is accurate, and easy to use. 
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APPENDIX 


NASTRAN COMPUTER INPUT DECK 
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TABLE 1. FREQUENCIES (CYCLES/SEC.) FOR CHANNEL VIBRATION 
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8. Beam along elastic axis v,0 

9. Beam along elastic axis v,0 

10. Beam along elastic axis v 

11. Beam along elastic axis v,0 


Natural Frequencies 


Warping 


Twist 


Bending 

1 

2 

3 

12 3. 

Included 


— 

’ 

— " 

— 

63.9 

255.5 

574.8 

No 

7.1 

14.2 

21.2 

— 

— 

— 

Yes 

15.5 

56.9 

125.8 

— 

— 

— 

Yes 

15.4 

56.8 

125.5 

67.1 

268.3 

603. 7 


No 

— 


— 

63.8 

254.4 

563.3 

No 

7.0 

13.3 

18.2 

63.8 

254.4 

563.3 

No 

— 

— 

— 

7.0 

12.9 

17.7 

No 

6.7 

12.7 

17.5 

66.6 

265.7 

588.1 

Yes* 

81.2 

134.7 

187.3 

59.0 

269.9 

591.0 

Yes 

— 

— 

— 

17.8 

63.9 

138.5 

Yes 

14.8 

54.1 

117.7 

66.7 

266.2 

589.1 


* Flanges are rigidly linked to the shear center. 









